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The 2D metal-insulator transition can be induced either by decreasing the carrier 
density or by increasing the spin polarization by applying a magnetic field parallel 
to the plane. Using experimental results for the shift in critical carrier density in 
an applied field, we show that the two-electron pair correlation function at short 
distance (A^r < 2) has a universal form along the transition line in the magnetic 
field-carrier density plane. This result provides direct evidence of the central role of 
Coulomb repulsion and the exchange-correlation hole in driving the metal-insulator 
transition. 
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The metal-insulator transition is a problem of great current interest 1 . In spite of intense 
investigation there is still no consensus on understanding the basic features of the problem 
and there is no microscopic theory. It is therefore useful to investigate phenomenological 
aspects of the metal-insulator transition. In this paper we point out an interesting empirical 
relationship between the dependence of the 2D pair correlation function g(r) on a magnetic 
field parallel to the plane of a two dimensional electron liquid and the freezing line in the 
phase diagram of the metal-insulator transition. 
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It is possible to regard the transition from the metallic state to an insulating state 
as a special case of freezing of the electron liquid. Phenomenological criteria in terms of 
pair correlation functions and static structure factors for the freezing of both classical and 
quantum liquids in the absence of impurities have been known for some time. The properties 
of the static structure factor S(k) near freezing show common features along the freezing 
line for a wide class of 3D interacting systems, and equivalent statements can be made in 
terms of the correlation function g(r) = 1 + J rf 3 k/ (2ir) 3 e* k,r {S'(k) — 1}. Hansen and Verlet 2 
established for systems of particles interacting through a Lennard- Jones potential in three 
dimensions that the Lindemann ratio along the melting line in the temperature-density 
plane is practically constant along the melting line. Hansen and Verlet also showed that 
on approaching the transition from the liquid state, the amplitude of the main peak in the 
static structure factor S(k) is very nearly constant ~ 2.85 ±0.1 along the crystallization 
line. This same value has been found for the peak in S'(k) along the crystallization line 
for both classical one component plasmas and systems of classical hard spheres 3 . Adler 
and Wainwright 4 obtained analogous results for two dimensional systems of hard discs. 
Generalizations of these ideas to quantum systems have been given, although the extra role 
of zero point motion makes the effect less clear-cut 5 . The essential point is that conditions 
for crystallization of the system as a function of external parameters correspond to critical 
values for short distance properties of the pair correlation function, or correspondingly the 
static structure factor. Consequently, insight into the freezing mechanism may be obtained 
by studying short distance properties of pair correlation functions as parameters are varied 
in the vicinity of the transition. This is clearly simpler than an explicit calculation of the 
equation of state, specifying full details of the system and its interactions. This can be useful 
for gaining phenomenological understanding in the absence of a microscopic theory. 

This characterization of the freezing line in terms of the strength of short distance cor- 
relations has been developed for pure systems for which the density distribution is uniform. 
The electron distribution in the systems exhibiting a metal-insulator transition is nonuni- 
form but since the transition is observed at low levels of impurities, the short distance g(r) 
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for kpr < 2 should be well approximated by the corresponding function for the pure system. 
Diffusion quantum Monte Carlo (DQMC) simulations 6 ' 7 give us the ground state properties 
of the 2D system without impurities including g(r). For the pure system the transition 
from liquid to solid occurs at extremely low densities corresponding to r s = 37 ± 5, where 
r s is the average electron spacing in effective Bohr radii. The preliminary quantum Monte 
Carlo calculations of Chui and Tanatar 8 which include the effect of charged impurities have 
suggested that impurities may stabilize a solid ground state at densities as high as r s = 10. 
Although a precise value of r s for the transition was not determined, this result does indicate 
that low levels of static defects in the form of a weak quasi-random potential may facilitate 
solidification. 

The metal-insulator transition occurs in very high mobility quasi-2D semiconductor de- 
vices with low carrier densities, r s > 5, (although very recently Altshuler et al 9 have dis- 
cussed data suggesting that the phenomenon may persist up to densities as high as r s ~ 2). 
At the very low densities of the Wigner crystal, r s > 37, the strong electron correlations 
generate a large area of near- zero electron density surrounding each electron 6 . This repulsive 
hard core region forms part of the electron's exchange-correlation hole, the complete density 
profile of which is known from the DQMC g(r). 6 The results in Ref. 6 establish that the 
radius of the exchange-correlation hole shrinks as the electron density is increased above the 
solidification density r s ~ 37. However, the hard core effect persists in the electron liquid 
up to a density corresponding to r s ~ 7. This is a factor some 25 times greater than the 
density at the Wigner transition. 

The DQMC numerical simulations in Refs. 6, 7 show that the exchange-correlation hole 
is stronger for spin polarized electrons than for unpolarized electrons. This is due to the 
additional exchange acting between the increased proportion of parallel spin electrons. In 
Fig. 1 we see that polarizing the spins at a fixed density corresponding to r s = 9 significantly 
expands the relative area of the hard core region in g(r). By r s = 7 the region of zero density 
around each electron has disappeared for the unpolarized system, but for the polarized 
system at the same density it is still present. 
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We have proposed 10 that the reported suppression of the metallic state by a parallel 
magnetic field 11 ' 12 , is directly associated with the expansion of the area of the hard core 
caused by the partial alignment of the spins in the magnetic field. We found that fields 
H\\ < 1 T destabilize the metallic phase. The enhancement in the exchange means that the 
critical impurity density n\ at the transition for the polarized system is smaller than the n\ 
for the unpolarized system at the same carrier density. 

Hamilton et al 11 reported that the metal-insulator transition boundary in p-GaAs is 
shifted by a parallel field H\ \ = 0.6 T from hole density p s = 7.5 xlO 10 cm' 2 top s = 12.4 xlO 10 
cm" 2 , or correspondingly, from r s = 9 to r s = 7. This is shown in Fig. 2. The transition 
boundary gives the critical magnetic field H c (r s ) as a function of carrier density. Moving 
along the boundary in the direction of increasing density, the value of H c increases, thus 
increasing the alignment of the electron spins. 

We connect magnetic field and spin polarization by estimating the degree of polarization 
as a function of a parallel magnetic field H\\ using DQMC data for the pure electron system 
taken from Rapisarda and Senatore 7 . In Fig. 3 we show the critical H c needed to fully 
polarize the electron spins, determined by equating the Zeeman energy gain to the energy 
difference, (g/j,B/h)H c (r s ) = [E p (r s ) — E u (r s )], where E p (r s ) and E u (r s ) are the energies per 
electron for the fully spin polarized and unpolarized states, respectively. We take (ga z ) = 1.1 
for GaAs 13 . For fields H\\ < H c (r s ) we use a linear relation to determine the degree of 
polarization. At r s = 7 a field H\\ = 0.6 T induces 50% spin polarization. 

We emphasize that the pair correlation function in the vicinity of the transition is not 
expected to be significantly affected at short distances by low levels of disorder. For this 
reason the DQMC g(r) for the defect free system can be used for the short ranged pair 
correlation function at the transition. Our main result is shown in Fig. 4 where we see for 
kpr < 2 that the g(r) on the transition boundary remains fixed all the way along the observed 
transition line. This indicates that for a given sample the transition is determined by a fixed 
functional form of the short-range <?(r), and is equivalent to a critical density profile of the 
exchange-correlation hole. The fixed functional form for g(r) = g(r, H c (r s )) suggests a fixed 

4 



functional form for S(k). This is reminiscent of the role of S(k) in characterizing the phase 
boundary in a wide range of liquid- "solid" transitions. It indicates that the metal-insulator 
transition is driven by electron-electron interactions and suggests that the transition is an 
intrinsic many-electron effect and not an artefact of material dependent properties. 

We conclude that the metal-insulator boundary for the sample as a function of spin 
polarization and carrier density is determined by a fixed functional form of the two-electron 
correlation function at short distances. By defining an exchange-correlation hole radius r c 
as the value of r at which g(r) exceeds a fixed value << 1 such as 0.01, then along the 
transition line as a function of parallel magnetic field H and carrier density p we have, 

This means that a knowledge of the dependence of g(r) on p and H permits us to predict 
the value of the gradient dH c /dp\ cr it-u n e of the phase boundary. 
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FIGURES 

FIG. 1. Two-electron correlation function g(r) taken from Ref. 6 at r s = 9 for unpolarized 
system (dotted line) and fully polarized system (dashed line), showing the effect of exchange 
enhancement. 

FIG. 2. Plot of phase boundary of metal-insulator transition as a function of parallel magnetic 
field and carrier density r s taken from Ref. 11. 



FIG. 3. Critical parallel magnetic field H c as a function of r s . 

FIG. 4. Correlation function g(r) along the observed metal-insulator transition boundary. Dot- 
ted line: r s = 9, unpolarized. Dash-dot line: r s = 8 for H\\ = 0.4 T. Dashed line: r s = 7 for 
H\\ = 0.6 T. The three curves are essentially identical, indicating a universality of the short range 
correlations at the transition. 
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